Realization of quantum walks with negligible decoherence in waveguide lattices 
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Quantum random walks are the quantum counterpart of classical random walks, and were re- 
cently studied in the context of quantum computation. A quantum random walker is subject to self 
interference, leading to a remarkably different behavior than that of classical random walks such as 
ballistic propagation or localization due to disorder. Physical implementations of quantum walks 
have only been made in very small scale systems severely limited by decoherence. Here we show 
that the propagation of photons in waveguide lattices, which have been studied extensively in recent 
years, are essentially an implementation of quantum walks. Since waveguide lattices are easily con- 
structed at large scales and display negligible decoherence, they can serve as an ideal and versatile 
experimental playground for the study of quantum walks and quantum algorithms. We experimen- 
tally observe quantum walks in large systems (~ 100 sites) and confirm quantum walks effects which 
were studied theoretically, including ballistic propagation, disorder and boundary related effects. 



In classical random walks, a particle starting from an 
initial site on a lattice randomly chooses a direction, and 
then moves to a neighboring site accordingly. This pro- 
cess is repeated until some chosen final time. This sim- 
ple random walk scheme is known to be described by a 
Gaussian probability distribution of the particle position, 
where the average absolute distance of the particle from 
the origin grows as the square root of time. First sug- 
gested by Feynman [H] the term quantum random walks 
was defined to describe the random walk behavior of a 
quantum particle. The coherent character of the quan- 
tum particle plays a major role in its dynamics, giving 
rise to markedly different behavior of quantum walks 
(QWs) compared with classical ones. For example, in 
periodic systems, the quantum particle propagates much 
faster than its classical counterpart, and its distance from 
the origin grows linearlywith time (ballistic propagation) 
rather then diffusively (2j. In disordered systems, the ex- 
pansion of the quantum mechanical wave-function can be 
exponentially suppressed even for infinitesimal amount of 
disorder, while such suppression does not occur in classi- 
cal random walks. 

In recentyears QWs have been extensively studied the- 
oretically [2J and have been used to devise new quantum 
computation algorithms 0]. Both discrete and continu- 
ous time QWs (DQWs;CQWs) flHB] have been studied. 
In DQWs the quantum particle hops between lattice sites 
in discrete time steps, while in CQW the probability am- 
plitude of the particle leaks continuously to neighboring 
sites. Both types of QWs have been studied theoretically. 
Experimentally, many methods have been suggested for 
the implementation of DQWs [see[2], but only a small 
scale system consisting of a few states was implemented, 
using linear optical elements _|7|]. For CQWs, a few sug- 
gestions have been made @, [J, yet only one experimen- 
tal method have been implemented by realizing a small 
scale cyclic system (4 states) using a nuclear magnetic 
resonance system[I3|. Such systems are difficult to scale 



to much larger configurations. Moreover, even at these 
very small scales, errors attributed to decoherence have 
been observed. 

Here we suggest a very different implementation of 
CQWs using optical waveguide lattices. These systems 
have been studied extensively in recent years but not 
in the context of QWs and quantum algorithms. We show 
that these systems can serve as a unique and robust tool 
for the study of CQWs. For this purpose we demonstrate 
three fundamental QW effects that have been theoreti- 
cally analyzed in the QW literature. These include bal- 
listic propagation in the largest system reported to date 
(~ 100 sites); the effects of disorder on QWs; and QWs 
with reflecting boundary conditions (related to Berry's 
"particle in a box" and quantum carpets [lJ,ll3|). Waveg- 
uide lattices can be easily realized with even larger scales 
than shown here (10 2 — 10 4 sites with current fabrication 
technologies), with practically no decoherence. The high 
level of engineering and control of these systems enable 
the study of a wide range of different parameters and ini- 
tial conditions. Specifically it allows the implementation 
and study of a large variety of CQWs and show experi- 
mental observations of their unique behavior. 

The CQW model was first suggested by Farhi and Gut- 
mann [6], where the intuition behind it comes from con- 
tinuous time classical Markov chains. In the classical 
random walk on a graph, a step can be described by a 
matrix M which transforms the probability distribution 
for the particle position over the graph nodes (sites). The 
entries of the matrix Mj t k give the probability to go from 
site j to site k in one step of the walk. The idea was to 
carry this construction over to the quantum case, where 
the Hamiltonian of the process is used as the generator 
matrix. The system is evolved using U(t) = exp(— iHt). 
If we start in some initial state \^in), evolve it under 
U for a time T and measure the positions of the result- 
ing state, we obtain a probability distribution over the 
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vertices of the graph. This is described by 

i-gf = -djlVj + 7j,j+iV>j+i + lj„j-tVj-u 



(1) 



where ipj is the wave function at site j, dj is the number 
of sites connected to site j (dj = 2 in the ID nearest 
neighbor case), and 7i,j(= 7j,i) is the probability per 
unit time for the transition between neighboring [14.ll5||. 
This mathematical formulation is effectively identical to 
the well known discrete Schrodinger equation used in 
the tight binding (Bloch ansatz) formalism in solid state 
physics It is used to describe the evolution of a wave- 
function on a periodic potential, which is essentially the 
propagation of a quantum particle on a lattice 3, 13] ■ 

An immediate implication for the correspondence be- 
tween QWs and these processes is that many of the ex- 
periments in solid state physics described by the tight- 
binding model could serve as implementations of QWs. 
However, such experiments deal with the macro-physics 
of the system and with overall observables such as con- 
ductance or transmission. Therefore, one can not mea- 
sure the specific spatial and temporal distribution of 
the electrons or photons wave-functions and the micro- 
physics of the system can not be directly observed. More- 
over, solid state systems contain many electrons which 
interact non-trivially and thus can not be described by 
the evolution equation of a single particle usually studied 
in QWs. Consequently, a qualitatively different experi- 
mental approach is needed in order to study QWs. Here 
we report such an approach using waveguide lattices. 

Recently, a new technique has been developed for the 
experimental investigation of periodic systems using op- 
tics. The salient feature of these experiments is that evo- 
lution of waves in time is also spread out in space, mak- 
ing it much easier to observe. This is achieved by using 
waveguide structures which are periodic on one dimen- 
sion (x-axis; see Fig. UK), but are homogeneous along the 
other (z-axis). In this way the wave propagation along 
the z-axis is free and corresponds to the evolution in time 
1 11 ] - Under appropriate conditions light is guided inside 
the waveguides and can coherently tunnel between them. 
The experimental setup and typical lattice parameters 
are described elsewhere (l8|. 

Light propagating in weakly coupled, single mode 
waveguides, can be described by [19]: 
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PjAj + Cj <j+1 A j+1 + Cjj-iAj-i (2) 



Here Aj is the wave amplitude at site j, f}j is the on- 
site eigenvalue, Cij is the coupling constant or tunneling 
rate between two adjacent sites i and j (for a periodic 
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Figure 1: (a) Schematic view of the optical waveguide lattice 
used in the experiments (see text), (b) Image of the output 
light distribution as recorded in the infrared camera, when 
the light is injected to a single lattice site at the input. 



main differences are that (1) the spatial modulation of 
the index of refraction in the x direction now plays the 
role of the tight binding potential, and the /3jS represent 
the propagation-constant eigenvalues of each waveguide 
in the lattice (2) the evolution at a given time can be 
observed by measuring the intensity distribution at the 
corresponding position in the z - axis [HI], since z = ct/n, 
where c/n is the speed of light in the medium. The ad- 
vantage of this system is the possibility to control the 
exact initial conditions for the light propagating inside 
the lattice. This is done by setting the width, the phase 
and the position across the lattice of the beam injected 
into the structure. In addition, this approach enables 
direct observation of the resulting wave-function by mea- 
suring the distribution of light intensity at the sample's 
output (fig. Wp>)- Furthermore, the temporal evolution 
of the wave-function can be observed by chan ging the 
sample length, or the initial conditions (e.g. 



One of the hallmarks of QWs on ordered lattices is 
their ballistic propagation [2j]. In order to observe this 
behavior, coherent light is injected into a single site in the 
lattice and the output intensity is measured. In Fig. [2]we 
compare the theoretical and the measured output distri- 
bution. The signature of ballistic propagation is clearly 
observed both at short and long propagation times (Fig. 
[2£i). Note that decoherence effects are negligible even 
after relatively long evolution in time, maintaining the 
detailed interference pattern predicted by theory (Fig. 
Gb)- Similar results, studied in a different context, were 
observed as early as in 1973 by Somekh et al. [2l| on 
small scales in structures similar to the ones described 
above. The propagating photons tunnel from the origin 
site to an adjacent site, and immediately start tunnel- 
ing to the next neighboring site. Through the tunneling 
between sites the photons accumulate a ir/2 phase, and 
an additional phase is accumulated continuously in each 



lattice Cij = C is constant), and z is the longitudinal lattice site j, at a rate given by (3j. The interference 

of all these waves depends on the phase accumulated in 
each possible path, and gives rise to the observed inten- 
sity distribution. This description is practically identical 
to the description of the QW, where the light intensity 



space coordinate. The description by Eq. [2]is completely 
analogous to the quantum description of non-interacting 
electrons in a solid crystal in the tight binding approx- 
imation, i.e. the discrete Schrodinger equation. The 
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Figure 2: (a) The theoretical prediction showing the ballis- 
tic evolution of the probability distribution of a CQW. The 
dashed lines correspond to the experimental measurements in 
(b). (b) The observed output pattern of light intensity after 
short (blue) and long (green) propagation in a periodic lat- 
tice. This well known pattern is one of the hallmarks of the 
ballistic propagation of QWs. (c) Output patterns of light 
intensity resulting from injection of light into two adjacent 
single waveguides (sites 42 and 43) of a disordered lattice. 
The different patterns observed demonstrate the high sensi- 
tivity of the QW to the initial conditions in this case. 



corresponds to the probability distribution of the quan- 
tum particle. Since the single photon and many photon 
problems are described by the same probability distribu- 
tion, experiments measuring light intensity are equivalent 
to performing a series of single photon experiments, from 
which the probability distribution is obtained. The prop- 
agation of more complex quantum states can be studied 
using correlated or entangled photons (see for example 
23]) . In this case the particle characteristics of the quan- 
tum walkers can be revealed by measuring two-photon 
correlation functions. 

When disordered lattices are used 2^, 24 1, a different 
behavior is observed. Accumulated random phases of the 
random walker lead to destructive interferences that in- 
crease with distance from the origin. As a result, after 
a short ballistic propagation, the tails of the distribution 
are exponentially suppressed leaving the probability dis- 
tribution exponentially localized to a small regime. This 



phenomena should be distinguished from a disordered re- 
lated decoherence. Decoherence is related to temporal 
disorder, which induces a loss of phase coherence and 
results in a transition into classical diffusion, character- 
ized by an expanding Gaussian probability distribution 
Spatial disorder such as used here leads to an 
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exponential (Anderson) localization [e.g. 
is a coherent interference effect . In the context of CQWs, 
such behavior was was found to be important for the ef- 
ficiency of quantum algorithms UJ, \23, |29j ■ 

QWs in disordered lattices are highly sensitive to the 
initial conditions. Fig. Eb shows two output patterns of 
light intensity resulting from the injection of light into a 
single waveguide of a lattice and similar injection to an 
adjacent site of the same lattice. The different patterns 
observed demonstrate the high sensitivity of the QW to 
the exact initial conditions. This serves as a unique sig- 
nature of the coherent nature of the QW, which is not 
present in the classical case. In addition these results 
demonstrate the effect of disorder on QWs, where in this 
case the disorder was introduced through randomizing 
the tunneling rate between sites (off-diagonal disorder). 
The tails of the distribution still show the ballistic com- 
ponent of the regular QW. However, additional strong 
peaks now appear near the origin. At later times these 
peaks evolve (on average) into an exponentially localized 
distribution, while the ballistic side lobes are suppressed 
(see [13] for detailed discussion). 

Several theoretical studies have been done on QWs 
with boundary conditions [3(1 l3l| , that give rise to com- 
plex self interference patterns. In Fig. [3] we show ex- 
perimental results of a QW with one reflecting boundary 
condition, compared with the theoretical analysis. A se- 
ries of measurements is shown (horizontal crossections) , 
where in each measurement light was injected closer to 
the boundary. The observed pattern results from the self 
interference of the incoming and reflected photons near 
the boundary, in agreement with theoretical predictions 
(30L I32T ] . Although these are limited observations show- 
ing results of a short time propagation, longer waveguide 
lattices could be used to study the more complex evo- 
lution at later times. For example, such behavior of a 
two boundary conditions system can be used for study- 
ing quantum carpets containing fractal patterns H, 23] ■ 

As an implementation of QWs, waveguide lattices carry 
some important advantages over other possible schemes. 
First, the technologies available for their fabrication or 
induction have reached a peak in recent years, enabling 
full control of every lattice parameters in ID and 2D 
[isL [331 ] , or limited yet real time control of lattice pa- 
rameters in 2D [34]. Second, waveguide lattices have ex- 
cellent structural stability, thus in practice decoherence 
due to noise is negligible. The optical wavelength in our 
experiments (using AlGaAs wafers) is around 1.5fim, the 
standard communication wavelength, and losses at these 
wavelengths are extremely small. This is highly impor- 
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Figure 3: (a) Measurements of the self interference patterns 
of QWs near a reflecting boundary. Horizontal crossections 
show the left half of the probability distribution of the QW, at 
decreasing input site position (vertical axis), where position 
marks the lattice left boundary, (b) Comparison to the 
theoretical analysis using the method of images [32]. 



tant for quantum computational tasks where coherency 
is essential. Third, effects arising from the interactions 
between different random walkers in other possible im- 
plementations are eliminated here, due to the bosonic, 
non interacting nature of photons. 

In recent years several quantum algorithms based on 
QWs have been suggested |35j. For realistic use of such 
algorithms one requires exponentially large systems. We 
note that as long as entanglement is not introduced, our 
system is limited to large but not exponentially large 
scale functionality. The lack of entanglement limits the 
number of the states of the system, which scales linearly 
with the number of waveguides. Our system, even with- 
out entanglement, can potentially implement QW algo- 
rithms, since quantum entanglement is not required for 
the algorithm implementation or its improved efficiency. 
Its only role in this case is to allow for a larger number of 
states (see, for example the discussion in [3a])- Some of 
the suggested QW algorithms have been shown to provide 
polynomial or even exponential speed up [37I . [38 ]. Un- 
fortunately, in all of the algorithms suggested so far the 
speed up of quantum over classical algorithms is achieved 
only when applied to high dimensional systems. Never- 
theless, our system can still be used to implement and 
study these algorithms in lower dimensions. 

In summary, we have demonstrated the strong corre- 
spondence between QWs and light propagation in waveg- 
uide lattices. This correspondence can be used to extend 
and interchange ideas and knowledge acquired in both 
fields (e.g. non-linear behavior [HJ in CQWs or entan- 
glement effects 

[sum 

in waveguide lattices). The high 
level of control, the accuracy, and the low decoherence 
rates achieved in waveguide lattices experiments provide 
a powerful tools for the study of QWs, and may prove 
useful in the implementation of QWs-based algorithms. 
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